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Abstract.  Aim  of  this  work  is  to  extend  the  results  obtained  in  a previous  study  on  the  magnetic 
confinement  and  stability  of  a quantum  degenerate  non-neutral  fermion  plasma.  This  extension 
consists  in  the  inclusion  in  the  previously  set  up  model  of  the  effects  of  the  exchange  forces,  and 
generalises  the  Thomas-Fermi  (TF)  approach  used  in  the  referenced  work  towards  a Thomas- 
Fermi-Dirac  (TFD)  statistical  description.  The  TF  model  has  not  only  been  used  extensively  and 
with  success  in  these  years  to  study  atomic,  nuclear  and  molecular  properties,  or  to  evaluate 
I'eatures  of  matter  in  extreme  conditions  such  as  low  temperatures  and/or  high  densities  typical 
of  astrophysics  and  inertial  confinement  fusion  experiments,  but  also  to  found  hydrodynamic 
theories  for  the  diffusion  and  stability  of  fermion  plasmas,  one  component  non-neutral 
degenerate  fluids,  plasmas  etc.  In  this  paper  an  equation  for  density  profiles  in  cylindrical 
symmetry  is  found,  from  the  semiclassical  kinetic  theory  of  quantum  gases,  which  takes  into 
account  the  effects  of  temperature,  average  velocity,  external  magnetic  field  and  quantum 
exchange.  Numerical  solutions  of  this  equation  for  the  case  of  complete  quantum  degeneracy  are 
given  and  comparisons  with  the  previous  results  are  carried  out. 


INTRODUCTION 

A many-body  charged  particle  system  in  which  there  is  not  overall  charge  neutrality 
can  have,  as  has  been  shown  since  the  early  ‘70s,  collective  oscillations  and 
excitations  as  well  as  shielding  and  screening  effects  [1];  this  many-body  collection  of 
charged  particles  is  therefore  termed  a non-neutral  plasma  (NNP).  In  certain  physical 
situations  the  NNP  is  in  such  conditions  that  also  global  effects  of  quantum  nature  are 
present  and  observable  [2,3,4].  These  conditions,  which  append  the  adjective 
“degenerate”  to  an  NNP,  are  realized  when  the  average  distance  between  neighboring 
particles  is  of  the  order  of  the  dimension  of  the  quantum  wave  field  of  the  particles 
themselves  [5]  (de  Broglie  length);  this  in  turn  depends  on  the  number  density  and 
temperature  of  the  plasma  itself:  the  lower  the  temperature  and  the  higher  the  density, 
the  higher  quantum  effects  are  and  more  and  more  degenerate  the  NNP  is.  Shielding  is 
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highly  affected  by  quantum  effects  and  the  screening  of  the  potentials,  with  its 
influence  on  the  spatial  dependence  of  density,  is  a main  subject  of  this  paper.  In  this 
work,  in  fact,  the  results  of  a previous  study  [6]  on  the  magnetic  confinement  and 
stability  of  a quantum  degenerate  fermion  NNP  in  the  presence  of  external  forces  are 
extended  by  including  the  role  of  exchange  potentials  on  the  particles;  the  extension 
consists  then  in  taking  the  description  from  a Thomas-Fermi  level  to  a Thomas-Fermi- 
Dirac  one.  Density  profiles  in  cylindrical  coordinates  and  symmetry  will  be  computed 
for  a completely  degenerate  electron  NNP.  Semiclassical  kinetic  theory  of  gases,  i.e. 
the  Boltzmann-Uehling-Uhlenbeck  (BUU)  equation,  and  quantum  statistics  in  the 
form  of  the  Fermi  statistics  are  the  starting  point  for  this  analysis. 


THE  MODEL 


For  the  case  of  an  electron  gas  the  BUU  equation  for  the  electron  distribution 
function  /can  be  written  (4,6|  as 

f - + V ■ ^ = |(/7>  (I  - >'/)('  - YA  )- #,  (I  - r/’)0  - Yf\  ))ga{f’,x)d^dx, 

d/  or  m dv  ^ 

(1) 

where  y = r.  Since  a stationary  solution  with  density  vanishing  at  infinity  is 

2m 

expected,  the  form  of  / is  known  to  be  a local  Fermi  distribution  [7],  for  which  the 
collision  term  of  eq.  (1)  is  vanishing;  it  is  more  convenient  then  to  make  use  of  the  H 
theorem  and  consequently  write  the  equation  in  the  structure  adequate  for  the 
logarithm  of  the  distribution  function  |7);  if  a finite  average  velocity  vo  niirst  also  be 
taken  into  account,  then  eq.  (I)  reduces  to: 
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where  c is  the  peculiar  velocity.  Introducing,  the  local  Fermi  distribution  function 
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into  eq.  (2)  the  following  equation  is  obtained: 
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If  eq.  (4)  is  to  be  satisfied  then  all  the  coefficients  of  successive  powers  of  the 
peculiar  velocity  between  square  brackets  must  identically  vanish;  hence: 


514 


(5a, b) 


mA  3v„ 


-cc:  ^ 
kT  3r 


= 0, 


dA 


v„  1 


dr 


I me" 

2kT^ 


3r 


ar  „ dA  A 

3r  9r  kT 


/ 

f 

\ 

F - w 

v„  — 

Va 

dr 

V 

) 

/ 

= 0 (5c, d) 


These  equations  are  as  general  as  eq.  (4)  but  are  indeed  extremely  useful;  in 
particular  eq.  (5c)  states  that  the  temperature  must  be  constant  and  eq.  (5d)  establishes 
an  equilibrium  relation  between  the  various  external  and  internal  forces  acting  over  the 
gas.  It  also  provides  a relation  between  the  chemical  potential  //  and  all  other  forces. 

We  consider  now  an  electron  gas,  in  condition  of  cylindrical  symmetry  with  respect 
to  z axis,  described  through  cylindrical  coordinates  r = [r,t5,z].  The  external  force  is 


given  by  a constant  and  uniform  magnetic  field  along  the  symmetry  axis,  B = Bz\  the 
remaining  internal  force  is  given  by  the  self-consistent  Vlasov  electric  field.  Under 
these  conditions  all  the  quantities  in  eq,  (5d)  will  depend  at  most  on  r alone.  This 
allows  a simple  form  of  vo  to  be  found  through  eq.  (5a);  taking  into  account  that  to 
have  stable  equilibrium  no  radial  velocity  must  be  present,  we  have  that 
v„  = v„.z -I- (0  X r , where  <d  = o)z  is  a constant  vector,  directed  along  the  symmetry 

axis,  representing  an  angular  velocity  for  the  gas.  Now,  inserting  this  expression  and 
the  Lorentz  force  into  eq.  (5d),  only  a scalar  equilibrium  equation  that  connects  the 
gradient  of  the  chemical  potential,  the  Lorentz  force  and  the  centrifugal  force  remains: 


kT — (ln^)-e(£,  +(orB)+  m^,co^r  = 0 
dr 

dV 


(6) 


Expressing  the  electric  field  as  — — ,withasb.c.  lim  , F(r)=  0,  we  get: 

dr  '■ 


/t(r,r)=  ju(^,T)+  eV{r)+  —mco{(0  - 0)^)r^  (7) 

where  /t(0,  f)  is  the  chemical  potential  for  an  uniform  and  homogeneous  electron  gas 
(4,5,8,9,10],  i.e.  the  value  of  the  chemical  potential  when  no  forces  are  present,  and 
eB 

co^  = - — is  the  cyclotron  frequency,  representative  of  the  external  magnetic  field.  The 

electric  potential  exchange  effects,  of  purely  quantum  nature,  are  taken  into  account, 
in  the  simplest  way,  through  the  Dirac  and  Slater  formalism  [8,9,10,1 1,12],  by  adding 
to  eq.  (7)  a suitable  term  proportional  to  the  cubic  root  of  the  local  particle  density 
(local  density  approximation,  LDA),  as  happens  in  the  usual  Thomas-Fermi-Dirac, 
TFD,  theory  of  atoms.  This  means  that  it  is  assumed  here  an  equal  number  of  upward- 
and  downward-  oriented  spins.  In  this  way  the  repulsive  effects  between  electrons  of 
parallel  spins,  generating  the  s.c.  “Fermi  hole”  or  “exchange  hole”  around  a generic 
test  particle,  are  considered;  however,  attractive  effects  due  to  other  kinds  of 
correlations  (f.i.  attraction  between  antiparallel  spins)  are  not  included.  Of  the  many 
relativistic  or  non-relativistic  LDA  models  [13]  the  non-relativistic  Dirac-Slater,  DS, 
one  has  been  chosen  because  of  its  simplicity,  both  from  a mathematical  point  of  view 
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and  from  the  physical  one.  The  fact  that  it  is  not  a relativistic  potential  is  coherent  with 
the  other  assumptions  of  the  present  model.  In  addition,  as  it  happens  with  other 
exchange  potentials,  the  DS  one  can  easily  be  related  to  variational  forms  and  inserted 
in  the  larger  theoretical  frame  of  density  functional  methods.  In  fact  the  DS  potential 
can  be  derived  from  the  following  density  functional: 


where,  in  SI  units,  C\. 


(8) 


IjL- 

4 Ake,, 


8;t 


. So,  the  DS  exchange  potential  is  simply: 


0/3 

y 


(9) 


It  must  be  observed  here  that  the  rhs  of  eq.  (9)  cannot  be  added  directly  to  eq.  (7), 
as  is  usually  done  in  the  TFD  method;  in  fact,  in  the  present  model,  eq.  (7)  is  naturally 
linked  to  the  condition,  different  from  the  usual  TFD  one  and  deriving  from  the  above 
stated  b.c.  for  V(r),  that  for  r-^0*  the  fermions  must  posses  the  characteristics  of  an 
uniform  and  homogeneous  gas;  because  n(f))*0  then  a very  important  normalizing 
term  must  also  be  added,  so  that  eq.  (7)  finally  becomes: 
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Inserting  eq.  (10)  into  the  Fermi  distribution  function  and  integrating  over  all  the 
velocity  space  via  Sommerfeld’s  Lemma  at  first  order,  one  gets  the  following  relation, 
generalized  to  take  into  account  finite  temperature  and  average  velocity: 
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where  /r„  = /r  ((),())=  is  the  s.c.  Fermi  energy. 

For  the  case  of  complete  degeneracy  (r=0)  one  gets  from  eq.  ( 1 1 ): 


an^'^  -k-hr^  +c  + dn'’^  =eV 


(12) 
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where  «(o)=r7„,  a--^ 
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I > 0 . For  r— >0"^  eq.  (12)  tends  to  an„^'^  + c = 0 which  is  the  usual 
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expression  for  an  uniform  and  homogeneous  electron  gas.  Also,  for  d—^0*  eq.  (12) 
gives  the  simple  TF  equation  in  absence  of  exchange  forces.  Requiring  now  the  self- 
consistency  of  electron  density  and  electric  potential,  imposing  that  eq.  (12)  satisfy 
Poisson’s  equation,  a non-linear  equation  for  only  the  electron  density  remains: 
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where  g=  — > 0 , with  particle  number  normalization  and  cylindrical  symmetry  b.c.: 


/i(f)=«(0)=«,„  lim^^,.  — = 0.  (14) 


ANALYTICAL  RESULTS 

A solution  vanishing  at  infinity  for  eq.  (13)  is  obtainable  if  and  only  if  [6]  its  rhs 
near  the  origin  is  negative;  this  means  that  to  have  confinement  it  must  be  gtj„  <46, 
independently  of  the  value  of  d.  This  important  result  shows  that  exchange  effects 
don’t  alter  at  all  the  s.c.  confinement  region  [6],  hence  they  don’t  alter  the  value  of  the 
confining  necessary  magnetic  field  too,  for  the  considered  gas  at  the  given  angular 
velocity.  Exchange  effects  do  indeed  alter  the  confinement  radius  for  the  gas,  in  the 
direction  of  extending  it  above  that  of  the  purely  TF  model,  as  is  to  be  expected,  as 
stated  above,  because  of  the  presence  of  “Fermi  holes’’. 


NUMERICAL  RESULTS,  COMPARISONS  AND  DISCUSSION 

Eq.  (13)  has  been  solved  numerically.  Figure  1 compares,  at  the  same  number 
of  confined  particles  per  unit  length,  magnetic  field  and  angular  velocity,  in  doubly- 
logarithmic  .scale,  the  TF  and  TFD  density  profiles.  The  role  of  the  repulsive  exchange 
holes  is  distributed  between  the  two  effects  of  lowering  the  density  at  the  origin  of  the 
TFD  solution  and  of  increasing  its  confinement  radius. 


CONCLUSIONS 

The  quantum  exchange  effects  between  parallel-spinned  electrons  have  been 
included  in  a semiclassical  model  of  a one-component,  non-neutral,  degenerate. 
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fermion  plasma.  The  plasma  rotates  around  its  simmetry  axis  in  a cylindrical 
configuration  and  is  magnetically  confined.  In  this  way  the  model  has  been  extended 


FIGURE  1.  Comparison  between  TF  and  TFD  models  for  N=270  in  ',  (O  = 4.2E+6  Hz  and  o)^  = 
8.0E+I4HZ. 

through  a local  density  approximation  from  a TF  description  level  to  a TFD  one,  as  it 
happens  f.i.  for  the  electronic  theory  of  atoms.  The  usual  expression  of  the  exchange 
potential  has  been  modified  for  the  special  needs  of  the  model.  The  equation  for  the 
density  profile  of  the  plasma  has  been  solved  numerically.  It  is  seen  that  the  s.c. 
confinement  region  is  not  altered  by  the  exchange  effects,  even  if  the  confinement 
radius  becomes  larger  because  of  the  repulsive  effects  of  the  “Fermi  holes”. 
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